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The theories of recoil-induced resonances (RIR) [J. Guo, P. R. Berman, B. Dubetsky and G. 
Grynberg, Phys. Rev. A 46, 1426 (1992)] and the collective atomic recoil laser (CARL) [ R. 
Bonifacio and L. De Salvo, Nucl. Instrum. Methods A 341, 360 (1994)] are compared. Both 
theories can be used to derive expressions for the gain experienced by a probe field interacting with 
an ensemble of two-level atoms that are simultaneously driven by a pump field. It is shown that the 
RIR and CARL formalisms are equivalent. Differences between the RIR and CARL arise because 
the theories are typically applied for different ranges of the parameters appearing in the theory. 
The RIR limit considered in this paper is qPo/Mcu q S> 1, while the CARL limit is qPo/Mui q < 1, 
£Nj , where q is the magnitude of the difference of the wave vectors of the pump and probe fields, Po is 

the width of the atomic momentum distribution and ui q is a recoil frequency. The probe gain for a 
probe-pump detuning equal to zero is analyzed in some detail, in order to understand how the gain 
arises in a system which, at first glance, might appear to have vanishing gain. Moreover, it is shown 
that the calculations, carried out in perturbation theory have a range of applicability beyond the 
recoil problem. Experimental possibilities for observing CARL are discussed. 
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I. INTRODUCTION 



O 

In recent years, there has been increased interest in spectral features that can attributed directly to the recoil atoms 
undergo on the absorption, emission or scattering of radiation. Among topics that have been discussed that fall into 
this category are Recoil- Induced Resonances (RIR) |l]-|| and the Collective Atomic Recoil Laser (CARL) [[To] [L3| . 
i' These processes appear to have much in common, although they are described quite differently. It is the purpose of 
this article to compare the RIR and CARL and to demonstrate that the RIR and CARL formalisms are equivalent. 
' The reason that this equivalence is not readily apparent is two-fold. Firstly, the RIR have been discussed using a 
density matrix approach in the Schrodinger picture, whereas CARL has been discussed using an operator approach 
(SJ ■ in the Heisenberg picture. Secondly, RIR and CARL have generally been examined for different ranges of the various 
parameters characterizing the atom-field interactions. In order to facilitate the discussion, it is useful to review briefly 
O ! the RIR and CARL. 

Since both the RIR and CARL relate to probe field absorption or gain in the presence of a pump field, it is perhaps 
best to recall the features of the probe absorption spectrum, neglecting all effects associated with recoil. Consider an 
ensemble of two-level atoms interacting with both a pump and probe field. The probe absorption or gain is monitored 
as a function the probe-pump detuning 6 for various pump field strengths, characterized by the pump field Rabi 
frequency xi ■ It is assumed that all collisional effects can be neglected and that the two-level atomic system is closed 
•^h | in the sense that the sum of ground and excited state populations is conserved for each velocity subclass of atoms. 
Moreover, it is assumed that the pump field detuning A from atomic resonance is much larger than any decay rates, 
Rabi frequencies, or Doppler shifts in the problem. In this limit, the probe absorption spectrum consists of three 
Q_i features, centered near S = 0, ±A Q. 

Of primary concern here is the structure centered near 5 = 0. In the absence of recoil, this line feature has a 
dispersion-like shape, is proportional to Xii and has a width of order of the excited state decay rate [jl5). Its origin can 
be traced to an interference effect involving processes in which two photons are emitted into previously unoccupied 
vacuum mode states . As such, it is linked to spontaneous emission, rather than a conservative exchange of energy 
between the pump and probe fields. The fact that the two-level system is "closed" plays a critical role here. Had 
the system been "open," the amplitude of this component would vary as Xi rather than x\ and its width could be 
determined by some effective ground state decay rate (e.g. inverse transit time) rather than the excited state decay 
rate 0. 

The RIR refer to a class of spectroscopic features in nonlinear spectroscopy that can be attributed to a recoil-induced 
"opening" of an otherwise closed, two-level system JjJ. In the presence of recoil, the atomic velocity is changed on the 
absorption or emission of radiation. As a consequence the sum of ground and excited state populations for specific 
velocity classes is no longer conserved. In the limit of large detuning A, the RIR can be interpreted as arising from 
Raman transitions between different center-of-mass states Q . If the pump and probe fields have propagation vectors 
ki and k2, respectively, then probe absorption occurs on a Raman transition between center-of-mass momentum 
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states |P) and |P— hq) and gain between states |P) and |P+?iq) , where q =ki— k2. Depending on the detuning S 
between the probe and pump fields, one of these processes is favored. For 5 < 0, there is probe gain, for 5 > 0, there 
is probe absorption, and for S — 0, the nonlinear probe absorption and gain vanish. In general, the RIR have been 
formulated as a stationary process in which the field achieves a steady-state value at each position R - gain occurs 
as the probe field propagates through the active medium. Steady state is achieved by assuming that there is some 
effective ground state decay rate or, alternatively, by assuming that the Doppler width associated with the Raman 
transitions, qPo/M, where M is an atomic mass and Po is the width of the momentum distribution, is larger than 
the recoil frequency, uj q = Tiq 2 /2M p,^|,^8|. The gain is normally expressed as dE^jdZ = a-E 2 , where i? 2 is the probe 
field's amplitude and a is a constant proportional to the atomic density. The point to note here is that the derivative 
of E2 depends locally on _E 2 . 

In contrast to the RIR, CARL has been formulated as a transient problem. It is assumed that the pump and probe 
fields are modes of an optical cavity in which the atoms are located. As a result of the atom-field interaction, the 
probe field builds up in the cavity as a function of time. In the limit of large A, one obtains an equation of the form 
dEijdt = j3 J Q f{t,t')E2{t')dt' , where (3 where is a constant proportional to the atomic density and / is a kernel. If 
qPo/M < ui q , the time derivative of E2 can depend on the past history of the field, rather than locally on the field. 
This leads to very different dynamics than those encountered in the RIR. One finds threshold conditions for probe 
gain that depend nonlinearly on the atomic density, as does the gain itself. CARL is a cooperative effect in the sense 
that gain occurs only for some minimum atomic density for most values of the detuning S. It is possible to have gain 
for 5 > 0, which is not possible for the RIR. On the other hand, if qPo/M ^S> uj q , the equation for dE^jdt becomes 
local and the gain characteristics of CARL become identical to those of the RIR. Thus, although both the RIR and 
CARL can be formulated for arbitrary ratios of qPq/Mui a , one can label the RIR limit as one in which qPo/M ^> uj q 
and the CARL limit as one in which qPo/M < uj q |lj|. This is the nomenclature that is adopted in this article. The 
Carl limit can be stated as, 2Po/frq < 1, which places it at or near the subrecoil limit of atom cooling. It should be 
noted that Verkerk and Bonifacio |n[] and Bonifacio et al. |l2j have shown that the CARL equations reproduce the 
RIR results in the limit that qPo /M ^> oj q ; however the approach they followed differs considerably from the one to 
be presented herein. 

It is important to recognize that it is possible for dE2/dt to depend locally i? 2 (t), even in the "CARL limit," 
qP /M < iv q . For example, if one includes ground state decay with rate T, a local dependence occurs if T ^> uj q ; as a 
consequence, in the presence of decay, the RIR limit can be stated as one in which qPo /M 3> w q or T uj q and the 
CARL limit as one in which qPo/M < uj q and V < uj q |I9|. Moreover, the build up of the probe field at very early 
times always varies linearly with the atomic density and can be viewed as a local process. In this sense, the theory 
of RIR is always valid at early times - the specific conditions are given in Sec. III. However, when q Pp /M < uj q and 
r<w 9 , the gain coefficient in CARL depends nonlinearly on the atomic density for times t > lu^ 1 pSJ^ , and the local 
dependence no longer holds. 

In Sec. II, the basic equations are derived and the RIR and CARL limits of these equations are obtained in Sec. 
III. An effective potential for ground state atoms is written which helps to categorize the RIR and CARL limits. The 
case of equal pump and probe frequencies is analyzed in Sec. IV using an effective five-level atom. When expressed 
in this form the calculation has a range of applicability beyond the recoil problem. The results are discussed in Sec. 
V. Special emphasis is placed on the distinction between "atom bunching" and "matter gratings." Experimental 
implications of the results are also explored. 



II. BASIC EQUATIONS 



The problem under consideration consists of a pump field and a probe field interacting with an ensemble of two 
level atoms . In some applications in CARL, it may be of interest to use a quantized description of the fields to follow 
the build-up of the probe field from noise, but, in the present discussion, the fields are taken to be classical. The 
pump field, 



Ei(R,i) = -ej ^i(R,i)e 4(kl - R " 0lt) +^ 1 (R,t)e" l(kl ' R -^ 1 * ) 



(1) 



has polarization e±, slowly- varying electric field amplitude i?i(R, i), propagation vector ki and frequency Q± = k\c 1 
while the probe field, 



E 2 (R,i) = ie 2 £; 2 (R,i)e l ( k2 - R -° 2 *) +£;*(R,t)e- l(k2 - R -^ 2 * ) 



(2) 



has polarization e 2 , slowly- varying electric field amplitude -E 2 (R, i), propagation vector k 2 and frequency £! 2 = fc 2 c. 
As a result of the nonlinear interaction with the fields, the probe field can be amplified. If the pump field is initially 



2 



much more intense than the probe field, as is assumed, pump depletion during the early stages of probe amplification 
can be neglected. Since the calculation in this paper is limited to the early stages of probe amplification, I take the 
pump field amplitude to be constant, 2?i(R, t) ~ E\. 

The dynamics of probe field amplification depends on the specific atom-field geometry. One can envision situations 
in which the probe field amplitude is constant in time but varies in space, or is constant in space but varies in time. 
The first case is the one considered generally in the MR, in which cw pump and probe fields interact with atoms 
in a cell or trap. The probe field amplitude increases in the direction of k 2 as it propagates through the medium, 
but is assumed to have evolved to a stationary state. The second case is the one considered generally in CARL, in 
which the fields correspond to field modes of an optical cavity and the atoms are located in the cavity. For a ring 
cavity, the probe field intensity is spatially isotropic, but increases in time Jjl|. In certain limits (to be noted below), 
the spatial gain coefficient of the RIR multiplied by the speed of light coincides with the temporal gain coefficient of 
CARL. In other limits, the dynamics of CARL amplification differs from the spatial build-up of the field in the RIR. 
To compare the RIR and CARL, it is convenient to adopt the cavity model and assume that the probe field amplitude 
is a function of time only, E 2 (R, t) ~ E 2 {t). All cavity losses are neglected, as is any ground state decay. 

Although CARL is referred to as a collective effect since conditions for CARL gain depend on the atomic density, 
each atom in the sample, on average, interacts with the fields in an identical manner. Thus, it is sufficient to write 
the Hamiltonian for a single atom interacting with the fields. The dependence of the field gain on atomic density is 
included automatically in the coupled Maxwell-Bloch equations. In the dipole and rotating wave approximations, the 
Hamiltonian for our system is 



P 2 Hco 
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H = 2M + T a * + H ^ M <} ei(VR_SV W + adjoint] , (3) 

where P is the atomic center-of-mass momentum, u> is the transition frequency between the ground state |1) and 
excited state |2) of the atom, a z = (|2) (2| - |1) (1|), <r+ = |2) (1|, 

*,(*) = -H^jp (4) 

is the Rabi frequency of field p, and p= (2\p\V) is a dipole moment matrix element. Terms related to spontaneous 
emission are not included in the Hamiltonian (||), for reasons to be discussed below. The Hamiltonian determines the 
time evolution of the atom, but the atomic evolution must be coupled to the field evolution via Maxwell's equations 
to arrive at a closed set of equations. 

The time evolution of the (complex) probe field amplitude is linked to the polarization of the medium which, in 
turn, is determined by the atom-field interaction. The medium's polarization can be expressed in terms of atomic 
density matrix elements as 

P(R,t) = N[pp 12 (H,t) + p* p 21 (R,t)}, (5) 

where N is the total number of atoms. As a result of the nonlinear atom- field interaction, the density matrix element 
p 2 i(R,t) can be written quite generally (see below) as 

/ 9 2 i(R,t) = ( 5 2 i(i;l)e l(kl - R ~ 0l * ) + ( 5 2 i(i;2)e l(k2 - R ~° 2 * ) + ^ p 21 (t; 1, n) e i(ki-R-o 1 t)+in(q-R+5i) ) (g) 

n/0,-1 



where 



q = k 1 -k 2 ; S = Q 2 -il 1 . (7) 



In the slowly varying amplitude and phase approximation, it follows from Maxwell's equations and Eqs. (^]) and (|^) 
that the evolution of the probe field is given by 

dE 2 (t) _ in a _ ca _ p(R t) e _ j(k2 . R _ n2t) (8a) 



dt €q 

_ iNn 2 



e 2 -p*p 21 (t;2), (8b) 



where P 2 (R, t) is the component of the polarization P 2 (R, t) varying as e^ k2 ' R n2t \ Combining Eq. (8b) with Eq. 
(i yields 
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d X2 (t) l Nd 2 n 2 

where 

d=\p-e 2 \. (10) 

An equation for p 2 \(t\ 2) is obtained by solving the Schrodinger equation for the Hamiltonian (|^). 
It is convenient to expand the wave function of the system as 

^(R,r,t) = ]T ^( R >*)^v( r ) (11) 

0=1,2 



where 



(12) 



uji = —co/2, uj 2 = cj/2, i?p = P 2 /2M, and Vv( r ) ^ s ^ ne atomic-state eigenfunction associated with state /i. Note 
that the momentum state amplitudes, A^(P,t), are defined in an interaction representation. The state amplitudes, 
a /J (P,i), in the "normal" representation are related to those in the interaction representation by 

a^P, t) = A P (P, t) e-^ 1 e~ lEpt ' h (13) 

and density matrix elements in momentum space are given by 

P^(P,P',t) = a M (P,t)a*(P / ,t) = A^V^AUV ,t)e-^it/* e -^ (14) 

where u!^ u — — lo u and -Epp/ = Ep — Epi . When the wave function (11^) is substituted into Schrodinger's 
equation with the Hamiltonian ([|), one finds that the state amplitudes A fi (P,t) evolve as 

ii(P,<) = -i X*„(t) exp[iA^i - iu k t - ik„ • Pt/M] A 2 (P+hk M ,t) (15a) 

0=1,2 

A 2 (P,t) = -i ^ XmW exp[-iA M t-iw fcM t + ik M .Pt/Af]i4 1 (P-ftk M ,t), (15b) 

0=1,2 



where 



is an atom-field detuning and 



A p = fi„ - w (16) 



u,, - hk 2 J2M (17) 



is a frequency associated with atomic recoil. In terms of these state amplitudes, the density matrix element p 2 i(R, t) = 
A 2 (R,t)Al(R,t) can be written as 

0u(R,t) = l - TJ1 f dPdP' 'A 2 (P \t)A\{P' \t) e -*-* e *{P-P')-R/* e -iE PP ,t/h (18) 

(27T?j) J 



Together, Eqs. (g), (|6j), (18) and ( p_5[ ) form a closed set of equations which can be solved numerically to obtain the 
time evolution of the probe field. 

To simplify the analysis of both the MR and CARL, it is assumed that 

|A 2 |«|A 1 |^|A|»fc Al u,| XAI |, 72 ; |5|«|A| (19) 



where 72 is the excited state decay rate and u is the most probable atomic speed. If inequalities (|19|) hold, both the 
RIR and CARL can be interpreted in terms of stimulated processes involving the pump and probe fields - spontaneous 
emission plays a negligible role. Both the pump and probe fields are turned on and brought adiabatically to their 
"initial" values at t = in a time long compared with |A| , but small compared with all other evolution times in 
the problem. Thus, at t = 0, the values for the Rabi frequencies are 
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Xi(Q)=Xi; Xa(0)=X2(0), (20) 
while the density matrix for the atoms is taken as 



g^(P, P'; 0) = a M (P, 0)a* (P', 0) = A^P, 0)A* (P', 0) = (27r?i) 3 y- 1 ^(P)(5(P - P') 8^8^, (21) 

where is the sample volume, W(P) is the initial momentum distribution, and the bar indicate an ensemble average 
pp| . This initial density matrix corresponds to a Wigner function W(R, P; 0) = W(P)/V, which is the product of the 
momentum distribution times a uniform spatial density (recall that this is a single particle density) . The calculation 
can still be carried out using an amplitude picture. The initial condition ( pl| ) is invoked once the density matrix 
element £> Mi ,(R, t) has been evaluated. 

When inequalities ( |l9| ) hold, it is possible to solve Eq.(15b) adiabatically for the upper state amplitude in terms of 
the lower state amplitude. Explicitly, one finds 

A 2 {P,t) = [XmW/ a 1 exp[-iA ll t-icj k t + ik IJ ,-Pt/M}A 1 {P-hk fi ,t), (22) 

where ojk = ^k t an d ter ms o f order |(Ai — A2) /Ai| and \(u>ki — UJ k 2 )/ U! k 1 \ have been ignored. Substituting this 
expression back into Eq. (|15aD, one obtains 



dA 1 (P,t)/dt = -i X ^X»® exp[iA vll t - icj k ^t - ik vil ■ pt/M] Ai(P+ftk V|1 , t), (23) 

v=\, 2 

where 

A„„ = A„ - A M ; = k„ - k M . (24) 

Since k„ M can take on the values [0,±q = ± (ki — k2)] only, it is clear that a solution to this equation can be written 
as 

00 

A 1 (P,t)= S n (P,t)A 1 (P-nnq,0). (25) 

n— — 00 

Combining Eqs. @, (|l|), @ and @, one finds that the density matrix element p2i(R,t) is given by 
p2i(R,t) = ^ 5^ XM ^ c ifk "' R ~ n " t V f "~"' )q ' R c~ u " 2 ~"' 2w 

fi— 1,2 n.n' — — 00 

x ydP^(P)e-^"- ,l ') p - qt / M S* n (P + nnq,i)S**(P + n^q,i). (26) 

Equation ( p6| ) proves that the general form for /92i(R, i) is correctly given by Eq. (^|). Extracting the coefficient of 
e z(k 2 R-n 2 t) p r0 vides us with the value of p2i(t; 2) needed in Eq. (^). Two types of terms in the summation vary as 
e j(k 2 -R-n 2 t) ^ ^ nose w ith p = 2 and n = n' and those with /i = 1 and n' = n + 1. Denoting the n = n' term as p~2i(i; a )> 
the n' = n + 1 as P2i{t\ b), and using the normalization condition |p(if 

1=53 /dR|A^(R,t)| 2 « /"dR |A!(R,<)| 2 = / dPW(P) |5„(P + riTiq, i)| 2 (27) 

one finds 



with 



M*; 2) a) + &1 (*;&), (28) 

fri(*5«) = ^ f] /dPW/(P)|S„(P + nftq,i)| 2 «^; (29a) 

n— — 00 

P2i(i;6) = ^- ^ e <(2n+iK*y dP^(P)e- lP ' q */ M 5„(P + n?iq,t)5: +1 [P+(n+l)^q,i]. (29b) 



n— — 00 



The first term represents linear dispersion of the medium. The second term, which is at the heart of the RIR and 
CARL, represents a scattering of the pump field off the atomic density distribution created by both the pump and 
probe fields. We now have all the ingredients necessary to derive the RIR and CARL equations. 
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III. CARL AND RIR 



In the large detuning limit, both the RIR and CARL equations are most easily derived using an effective Hamiltonian 
for ground state atoms. An effective Hamiltonian of the form 



Heff = ^ + ^(\xi\ 2 + \X2(t)\ 2 ) + ± [xixUt) e^ R+5 *> + xMt) e-^ R ^} (30) 



leads to the evolution equation (|23|) for the ground state amplitude Ai(P, t). The second term in the Hamiltonian is 
the spatially homogeneous light shift potential and is not of interest here. The third term is the potential formed by 
the pump and probe fields which gives rise to RIR and CARL. The depth of the potential, U = fiXiX^^)/ is a 
critical parameter. If the energy spread of the atoms 5£ is much larger than U, any "bunching" or spatial modulation 
of the atomic density will be negligibly small, independent of the detuning 5. This is generally the case at thermal 
energies, making it all but impossible to observe RIR and CARL using atoms in cells at or above the Doppler limit 
of laser cooling (see discussion in Sec. V). Assuming that 5£ < U, there are still two limiting cases to consider. If 88 
is much greater than the level spacings in the potential U, then the motion of the atoms can be treated classically. 
By approximating the potential as that of a harmonic oscillator and assuming that the spread AP is of order hq, 
one can estimate that the validity condition for the motion to be treated classically is \xiX2/ Acu q \ ^> 1. This is the 
limit considered by Bonifacio and coworkers and will be referred to as CARLSC, with the "SC" denoting the 
semiclassical limit. A theory of CARL in which the atomic motion is fully quantized has been given recently by Moore 
and Meystre It will be referred to as CARLQ, when it is necessary to distinguish between CARLSC and CARLQ. 
The theory of RIR has been developed using a quantized description of the center-of-mass motion; consequently, the 
RIR and CARLQ theories should produce identical results, and the RIR and CARLSC theories should produce the 
same results when |xiX2 / | 3> 1. It should be noted that, although the atomic motion in the optical potential is 
treated classically in CARLSC, the gain in CARLSC results from recoil-induced processes. This point is discussed in 
more detail in Sec. V. 

In calculating expectation values of operators that are diagonal in the internal states, one can neglect contributions 
from the excited state since the excited state population is assumed to be negligibly small. The CARL equations are 
derived using the Heisenberg equations of motion while the RIR are derived using a density matrix approach. I now 
show that these methods yield equivalent results, as they must. 

A. CARL 

To make connection wit h CA RL, one defines a Heisenberg operator, "bunching parameter," B(i) = e 1 ^' 11 ^ . Using 
Eqs. (U), ©, ©, and (p9b|), one can write the average value of the bunching parameter as [M 

(B(t)) = (e l * R ^= J dPdP' p n (P,P';t) (P'|e' iq,R |P) (31a) 
= 1 J dXLdPdP' A 1 (P,t)A* 1 {P', t) e iq-R e <(P-P')-R/R e -iE PP ,t/n ( 31b ) 
= J dP A 1 (P,t)A* 1 (P + hq,t) e -iEp, ]P+ n^t/h ( 31c ) 

= / dP S n (P,t)S* l ,(P' + Uq,t) Pll [P -nhq,P' - (n' - l)fiq; 0] e _iBj '-l p + R -i t/ ' i (31d) 

n,n f = — oo 

oo „ 

^V- 1 Y e i(2n+1)u, " t / dPW(P)e" lP - rlt/M S n {P + nn( l ,t)S; i+1 [P + (n+l)h( l ,t} (31e) 

71 — — OO 

AVe i5t 
Xi 

where the brackets denote a quantum-mechanical average and the bar an ensemble average. Therefore, by combining 



Eqs. @, @, @, and @, one finds 

d X 2(t) _ nd 2 n 2 f X2(t) , Xie l5t 



-i- 



dt 2heo 



Bit)) , (32) 



G 



where n — N/V is the atomic density. This is one of the basic CARL equations in the limit of large A E3]. The 
remaining CARL equations, obtained from the Heisenberg equations of motion with the Hamiltonian (30) are 



dB/dt 



q-P(*) 

M 



u q B(t); 



dP/dt = -A XlX *( t )e iSt B(t) ~ X2 {t) X {e~ m B\t)]. 



(33a) 
(33b) 



These equations do not form a closed set since the equation for (dB(t)/dt) involves (P(t)B(t)}. One must generate a 
series of equations for these higher moments or use some approximation techniques (such as assuming classical motion 
in the effective potential) to obtain a solution to the equations. The linear dispersion parameter a = ^u/a* resu lts in 
a shift of the probe field frequency. By setting X 2{t) — X2{t)e l8t in Eqs. ( |32| , [33"| ) , it is easy to see that the detuning 5 
enters only in the combination 8 + a. Thus, one can redefine the detuning to include the dispersion shift and replace 
Eq. © by 



dX2(t) 
dt 



. nd 2 n 2X ie M 

1 2he A {B{t)) 



(34) 



To examine the small signal gain regime, one can solve Eqs. (|33|) in perturbation theory. To zeroth order in the 
fields, 



P(t) = P(0) = P; 



B(°\t) =exp<ji 



B<°>(t) 



— exp < —i 



q P 

M 

./qP 



uj q t\B(0) =B(0)exp\i 



q P 

M 



+ UJ q )t 



(35a) 
(35b) 

(35c) 

_ p iq R(0) — „iq R 



where the last equalities in Eqs. ( p5qj35c[ ) follow from the commutation properties of P and B(0 ) 
Note that P = P(0) and R = R(0) are normal Schrodinger operators. Substituting Eqs. ( 35b ,35cj) into (|33b|) , one 
finds, to first order in IX1X2I7 



M 



t}B\0) = S t (0)exp 



q p 



pW(t) 



dt' 



-X2{t')x*i e- l (^ +s+ ^y St (0) 



(36) 



and, when this result is substituted in Eq. (33a), one obtains 



B«(i) = 



2u> q e~ 


-iSt 


A 








*/' 


dt" 


JO 




2uj q 






f dt 


A" J 









*/' 


dt" 


Jo 





dt 1 e^ l{ ^ +S -^ (t - tl] 

-Xa(*")Xi e-i^+^y Bt(o) e < s kr+' 5 -<^>' B(0) 



' I XiX^(t")exp i; (^ +5 ) ( *' +t " ) e^( t '- t ")[ J B(0)] i 
-X2(t")x* exp l (^ +,5+ ^)( t '-*") 



(37a) 



(37b) 



where the commutation properties of B(0) [or B'(0)] and P and the relationship B'(0)B(0) = 1 have been used 
to go from (37a) to ( |37b| ). The bunchin g par ameter B(t) = B^(t) + B^(t) can now be aver aged with the initial 
density matrix (2l|). The average of EqJUb)) for B^(t) vanishes as does the first term in Eq. $Wj) for B^(t). On 
averaging the remaining term in Eq.(37b) for B^(t) with the initial density matrix ((2^), interchanging the order of 
integration, and carrying out the integration over i',one obtains 



(B(t)) =- 2X \ ~ J dPV^(P)^dt'x 2 (Oexp ,; (^ +5 ) (t -*' ) sinK(t-<')] 



(38) 



I will return to this equation after showing that an identical equation is reached using a density matrix approach. 
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B. RIR 



The RIR are usually calculated in the context of a density matrix approach. From Eq. (31a), it follows that 

(B(t)) = J Pll (P,P + ?iq;<)dP, (39) 
which gives the alternative form for Eq. (|34j), 

■Jpu(P,P + tuKt)dP. (40) 



dx2(t) _ nd 2 Q 2 Xi e 
dt ~ % 2he A 



To complete the RIR equations, one uses the Hamiltonian (|3C|), along with Eqs. (p])-(|l4|), to obtain density matrix 
equations 

0p u (P, P'; t)/dt = -iu> PP ,pn (P, P'; t) 
■ X*iX2(t)e- lS! 



A 

■xix*2(ty st 



- [p n (P + hq, P'; i) - p n (P, P'-hq; t)] 
[ Pll (P ~ hq, P'; t) - Pll (P, P'+hq; t)] , (41) 



where ojppi = Epp> /%. It is easily verified that equations for the quantities d (B(t)) /dt = j pn(P, P + hq;t) dP and 
d(P(t)) /dt = JPpn(P,P;t) dP, obtained using Eq. (^l|), are identical to Eqs. (p3|). Thus, the RIR density matrix 
equations are totally equivalent to the corresponding operator CARL equations. 

To evaluate (B(t)) in the perturbative limit using Eq. (j39|), one sets P' = P + hq in Eq. pi] ) and replaces the 
density matrix elements in the right hand side of that equation by their zeroth order values, 

pff (P, P'; 0) = (2Trh) 3 V- 1 W(P)S{P - P'), (42) 

to obtain p0|| 

dp[\ ) (P,P + TU%t)/dt = -»Wp|P + n«i|/»n(P,P + &q;i) 

_- XiX2(*)e ' [M/(p + hq) _ w ^ (43) 

This equation is consistent with the RIR picture of Raman transitions between center-of-mass momentum states 
differing by hq. Integrating Eq. ( (43| ) over P and t, and using Eq. ([39]), one reproduces Eq. ([58]) for (B(t)). Since 
(B(t)} ~ as w ? ~ 0, the probe gain in the RIR and CARL is a recoil-induced effect. 
The value of (B(t)) given by Eq. (gg), which determines the small signal gain, depends critically on the ratio |lj| 

- - qP ° 2F ° (44) 



M LO q hq 

where Pq = Mu is the width of the momentum distribution and u is the most probable atomic speed. For r 1 
("RIR limit"), the integrand in Eq. ( |38| ) is rapidly oscillating except when t' f=a t, allowing one to approximate the 
integral as 

{m) ^ _ 2u qX le-™ X2 (t) J dp w{p) J* dt , {t _ t , )e< ^p +5)(t _ (45a) 

^^^ j/ dPW(P) £ dt>e<^^-t'A /d6 (45b) 

= i ^ipm (^)d{[ dl Wq (P 0X ) J'" dye<^)) /d8, (45c) 

where it has been assumed that the momentum distribution can be written as the product of a one-dimensional, 
symmetric distribution, W q (P q ), and a two-dimensional distribution, W±(P±), for momenta Pj^ transverse to the q 
direction. Equation (45) is the general RIR result, in which the bunching parameter depends locally on the field Xiip)- 
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By combining Eqs. ( [45| ) and (|34|), one sees that, in general, the exponential build up of the field is not linear with 
time. However, if qut 3> 1 and 6/qu < 1, the expression for (B(t)) reduces to 



(£(*)> = 



2nuj qX le- M X2 {t) M 



xd 



A V 1 

dW q (P q = gM/g) 
d<5 



dxW q (P x)V [x+—\ } /dS, 



qu 



(46) 



where V indicates a principal value, enabling one to combine Eqs. (ktq) and (34) to obtain the small signal gain [III 



nd 2 n 2 irq\ X i\ 2 dW q (6M/g) 
2e A 2 dS 



(47) 



The gain depends linearly on the atomic density and the pump field intensity. The gain is positive for S < 0, negative 
(absorption) for 6 > 0, and vanishes at 6 = 0. This is a "single particle" gain in that each atom contributes separately 
to the gain and there is no threshold condition for gain that depends on atomic density. 



The situation changes in the CARL limit r -C 1. In that limit one can replace W(P) in Eq. (38) by <5(P) to obtain 



(#(*)> = - 



2*i e 



-iSt 



dt' X2 {t') exp*(*- t '> sin [cj,(t - *')] 



(48) 



The bunching parameter now depends nonlocally on the field amplitude X2(^')j that is, it depends on the past history 
of X2 it) ■ Equation ( 48 ) is equivalent to the differential equation 



(B{t)) /dt 2 = -u? q (B(t)) 



-iSL 



(49) 



subject to the initial conditions (-B(O)) = d(B(Q)) /dt = 0. The small signal behavior is determined by the coupled 
equations (49) and (|34|). These equations have already been analyzed by Bonifacio and coworkers for CARLSC with 
neglect of the to 2 term in Eq. ([l9|) jl0| and by Moore and Meystre [jl3| including this term. Gain occurs if one of the 
roots of the cubic indicial equation, obtained from the coupled equations, ( f49| ) and (|34|), has a positive real value. In 
terms of the quantity 



Q 



nd 2 n 2 |xi|' 
hA 2 e a uj 2 



the indicial equation is 

and the condition for gain is [jjj 



+ iSs 2 + lu 2 s - ik)\{u q Q - 5) = 0, 



9 (6/u q ) - (S/w q y + 3 + (S/uj q y <27Q/2. 



3/2 



(50) 



(51) 



(52) 



If Q <C 1, gain occurs in the range \S/uj q + 1| < \/2Q and the gain coefficient equals y 2Q — (5/uj q + 1) u q /2. On 
the other hand, for Q > 1, gain occurs for 5/uo q > 0, in contrast to the RIR limit. Moreover in both cases, the 
gain depends nonlinear ly on the atomic density and there is a threshold condition for all values of S/uj q ^ —1. As 
such, CARL is a collective effect in the sense that gain does not occur for atomic densities below a certain critical 
value. For 5/uj q = —1, which is the resonance condition for Raman transitions between center-of-mass states having 
momenta P = and P = hq, there is gain irrespective of the value of Q. For Q> 1 and < 1, the gain varies as 

Q^I^LOq cos(7r/6) = (\/~3/2)Q 1 / 3 uj q , and the probe field undergoes exponential gain linear in time for t > (Q 1 / 3 ^) -1 . 



(53) 



For early times, one can evaluate X2(^') in Eq. (]3q) at t 1 — 0, and combine Eq. (|38|) with Eq. (M) to obtain 



l + ti 



dPW(P) / dt / drexp 



<^+ 5 ) T sin( Wg r)|x2(0). 



As long as the magnitude of the second term is much less than unity, the build-up of the probe field depends linearly 
on the density and the dependence of dxi/dt on X2 is approximately local. In this sense, the RIR limit is always valid 
for sufficiently small t or atomic density (since Q is proportional to the density) . It follows from Eq. (^3|) that the RIR 
limit is always valid if Quj 2 t 2 /2 < 1. If W(P) - S(P) and \S\ t,u; q t < 1, the RIR limit is valid when Quj^t 3 /6 < 1. 
The condition Qco 3 t 3 /6 <C 1 is less restrictive than Qu> 2 t 2 /2 < 1 for Q > 1. 







IV. DEGENERATE PUMP AND PROBE FIELDS, 6 = 



An interesting limiting case is one in which W^(P) ~ S[P) (CARL limit), and the pump and probe field frequencies 
are degenerate, fii ss £1 2 = f2; 6 = f2 2 ~ f2i = [recall that the detuning has been redefined to include the shift 
arising from linear dispersion]. When 5 = 0, the probe gain vanishes identically in the RIR limit r » 1, but grows 
exponentially for CARLSC. For CARLQ, it follows from Eq. J52] ) that the probe gain is exponential provided that 

Q > ITS (54) 

This qualitative difference between the RIR and CARL limits is reason enough to consider the S = case in some 
detail, but it is not the only reason. The existence of exponential gain when 5 — is surprising at first glance. It would 
seem that processes in which a pump photon is absorbed and a probe photon emitted would be exactly cancelled by 
processes in which a probe photon is absorbed and a pump photon emitted, owing to the symmetry of the interaction 
when W(P) ~ S(P). It is the purpose of this section to investigate the origin of exponential gain in CARL when 
5 = 0. Calculations are carried out in a perturbative limit, that is, to lowest order in the pump field intensity. As is 
shown below, the results are also applicable to a wider range of problems. 

In perturbation theory, there are five states that enter the calculation, starting from atoms in their ground state 
having P = 0. The relevant states are \g;P = 0) , |e;P = Tiki) , |e;P = Tik 2 ) , \g;p = ±27iq), having energies 0, 
h(ui + u>k), + u>k)i an d TiWq, respectively (recall that u>k = tOki ~ L °k 2 )- It is convenient to relabel these states as 

\g; P = 0) = |0) , |e; P = Tiki) = |1> , |e;P = ftk 2 ) = |-1) , \g;p= ±2hq) = |±2) . (55) 

The energy levels associated with these states are shown in Fig. [|. The pump field drives the |0) to |1) and | — 1) 
to |— 2) transitions while the probe field drives the |0) to | — 1) and |1) to |2) transitions. The Hamiltonian for the 
system can be obtained by expanding the Hamiltonian (^) in a momentum-state basis for the subspace (^5j). Using 
the relationship (P| e ik R |P') = 5(P - P' - hk), one finds 

H = h(Lu + Wfc )(|l) (1| + hi) (-11) + hu q (\2) (2| + |-2) (-2|) 
+Ti {xi e- int (|l) (0| + 1-1) (-2|) + adjoint} 

+h { X2 {t) e~ mt (\l) (2| + 1-1) (0|) + adjoint} (56) 

Note that this level scheme could equally well describe an atom in a Stark field, driven by circularly polarized pump 
and probe fields. The energy Tiuj q would then correspond to the Stark shifts of the m = ±2 ground state, Zeeman 
sublevels. 

Since decay is neglected, the calculation is most conveniently carried out using state amplitudes rather than density 
matrix elements. No ensemble average is needed here since we start in an eigenstate of momentum, P = 0. It is 
convenient to work in a field interaction representation in which the state amplitudes <2j(P, t) (j = —2, 2) are written 
as 



a (P,i) = 
oi(P,t) = 
o_i(P,*) = 
a±2(P,*) = 



with density matrix elements given by 



[(27rh) 3 /V}^ 2 ao(t)S(P) 
[{2nh) 3 /V} 1/2 ai {t) e- int 6(P - Tiki) 
[(27rTi) 3 /U] 1/2 5_i(t) e- int 5(P - Tik 2 ) 
[{2i,hf/V] 1 ' 2 a ±2 {t)5{PTK<il 



(57a) 
(57b) 
(57c) 
(57d) 



Pn,v{t) = a^{t)a* v {t) (58a) 

P**(R, t) = J dPdP'a^P, t)a* u (P', t) e '( p - p ')' R / fi . (58b) 

In Eq . (pal) for the probe field evolution, one need the component of polarization P 2 varying as e l ( k 2-R--^t) From 
Eqs. @and @, it follows that 

P 2 (R, t) = n[po,-iP-i,o(t) + P2,m, 2 (t)} e t(k2 ' R - m) (59) 
where the p's are dipole matrix elements. For the level scheme corresponding to the states (|55|), 
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P0,-1 = P2,l = Peg = P*- (60) 

Thus, the problem reduces to calculating p_i o(i) and pi,2(t) to order |xi| 2 X2- 



-1> - { |1> 




|0> 

FIG. 1. In perturbation theory, the relevant momentum states can be represented as an equivalent five-level system, interact- 
ing with two fields as shown. The initial momentum distribution is taken to be a delta function and the detuning is S = 0. State 
|0) corresponds to |gr;P = 0) , state |1) to |e;P = ?iki}, state — 1) to |e;P = 7ik2) , and states |±2) to |g;P = ±2frq). Although 
derived for the recoil problem, the conclusions reached in the text for this level scheme are applicable to a any problem where 
a similar level scheme is encountered. 

Before undertaking this calculation, it is useful to obtain an expression for the time evolution of the probe field's 
energy density W 2 =eo \E 2 \ 2 /2. Using Eqs. (§a|), (|), @, and @), one finds 

dW 2 /dt = innp*-e 2 EZ(t)(p- lt a + pi, 2 )/2 + c.c. (61a) 
= -inhflxl(t)(p-ifi + p\a) + c - c - (61b) 
= nhVL [p2.i - (p-i,-i + P-2,-2)] ■ (61c) 

The last line follows from the density matrix equations of motion for the Hamiltonian (^6|) and can be given an 
obvious physical interpretation. Population of state |2) implies gain on the probe field, while population in states 
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| — 1) or |— 2) implies loss for the probe field. The time rate of change in energy density is simply the difference 
[f>2,2 — (p-i,-i + P-2,-2)] multiplied by the product of energy, Ml, gained or lost by the probe field in each elementary 
process and the atomic density. In the adiabatic limit considered in this paper, the excited state population is negligible 
and one has 

dW 2 /dt ~ nM}(p 2 , 2 - P-2,-2). (62) 

One might expect that p 2 ^ 2 = p— 2,-21 owing to the symmetry of the level scheme, but we will see that this is not the 
case. 

Assuming that the detuning of the fields A from the ground to excited state atomic resonance is sufficiently large 
to adiabatically eliminate the excited states, one uses the Hamiltonian (|5^) to show that the excited state amplitudes 
are given by 



a iW ^-ooW 

a-iW 7T- a oW 



X2(t)~ M 

—£-a 2 {ty, 

Xi(t)~ ... 



(63a) 
(63b) 



and that the state amplitudes do(t), d± 2 {t) evolve as 



da (t)/dt = -ix*ai(t) - ixz(t)a,-i(t) 
da 2 (t)/dt = -iuj q d 2 {t) - iX2(*)5i(i); 
dd- 2 (t)/dt = —icj q d-2(t) ~ ix*5_i(t). 



(64a) 
(64b) 
(64c) 



It is a straightforward exercise to solve Eqs. ( p4[ ) in perturbation theory to third order in the fields, starting from 
o,o(0) = 1. One finds 



4 0) w 

a«(*) 

4 2) (i) 



1 

Xi 
A 

Xa(t) 



—00 
t 



—00 
t 



A 

dt l XlX2(t') iu a ( t -t') 

A 



% I dt 

' —00 



,IXi| 2 + | X2 (t')| 2 



A 



(*) 



a?l(*)~^4 9) (*) + fa«(t), 



(65a) 
(65b) 

(65c) 
(65d) 
(65e) 
(65f) 

(65g) 
(65h) 



where the superscripts denote the order of the fields. 

Consider, first, Eq. (|62|) for the probe field intensity, which depends on 



d(p2,2 ~ P-2,-2) /dt 



af(t) - a^(t) )/dt 



(66) 



By inspecting Eqs. (65d 65c), one can understand the manner in which 



a 2 2 \t) 



can grow more rapidly than 



~(2) 



(t) 



Suppose X2(t) acquires a positive, time-dependent phase as a result of the atom-field interaction. In this case, the 



quantity X 2 (t')e lu ' qt app eari ng in the integrand of Eq. ( |65d[ ) varies more slowly than the quantity X2(t')e lbJqt appearing 
in the integrand of Eq. ( |65e| ). As a result, state |2) population builds up more rapidly than that of state |— 2), leading 
to probe gain. In other words, the nonlinear phase modulation of the probe field effectively favors the |0) to |2) 
transition over the |0) to |— 2) if the phase is positive. 
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To examine the phase and amplitude build-up of the field, one can combine Eqs. (fa|), §), and © to 

obtain 



dX2 
dt 



ifld 2 

1 

2e h 

ifld 2 

1 

2e h 



Pl,2 







ft 



4 1} (4 2) ) ! 



^'x2(y'M) sin [H i -y')]- (67) 

Of course, the integral equation (|67j) is equivalent to the coupled equations (|34|,[58|) with S — 0. However, this form of 
the equations is convenient for obtaining the early time development of the field and for determining the conditions 
under which the field undergoes exponential gain. By solving Eq. ([37]) iteratively, one finds 



X2(t)/Xa(0) = [1 + iQu q t - {Quo q t) 2 /2] - iQsm(u q t) 



- (Q 2 /2)[4cos(w q t) + uj q t sin(w g t) 



o[Q 3 Kt) 2 )]. 



(68) 



For uj q t <C 1, X2 ~ 1 + i{Quj q t) 3 /6 



,t J /6 ^ ear iy times, the phase is positive, favoring probe gain. Whether or 



not the field continues to grow depends on the value of Q. If Q <C 1, then X2 is slowly varying compared with sin(uj q t) 
and the integral in (pTn can be evaluated asymptotically to yield 



which implies that 



dX2 
dt 



X2(t) 



iuj q Q[l - cos(w 9 i)]x2, 



JQ[uj q t— sin(oj,t)] 



X2(0) 



(69) 



(70) 



For Q <C 1, the intensity of the probe field remains approximately constant. In some sense, this can be viewed as the 
MR limit, since the time development of the field is approximately local, and there is neither absorption nor gain. 
On the other hand, for Q > 1, the field build-up occurs sufficiently rapidly to insure that 



(ha - ^-2,-2) = -2Re[X2(*)(P-i,o + Pi,2)} 



4 Re 



\xi\ 2 xm 

A 2 



dy'x 2 {t')sm(t-t') 



(71) 



remains positive for all t. In this case, there is exponential gain for the probe for times ujqQ^l^t > 1. 

In summary, the probe gain that occurs for 5 = and large A is clearly not a single particle effect. It is more 
closely related to a "propagation" effect in which the phase modulation of the probe field produced by the nonlinear 
atom-field interaction drives the probe gain. 



V. DISCUSSION 



It has been shown that the density matrix-RIR and Heisenberg operator-CARL formalisms lead to equivalent 
equations. The RIR and CARL limits refer simply to different regions of parameter space of these equations. For a 
given experimental situation, one must determine whether one is in the RIR limit, the CARL limit, or neither limit 
(as is most often the case). The experimental implications of the RIR and CARL are discussed below, but first I 
would like to discuss the distinction between the terms " matter grating" used in discussions of the RIR and " atomic 
bunching" used in discussions of CARL. 

The term "matter grating" refers to spatially modulated atom distributions resulting from a nonlinear atom- field 
interaction. The term "atomic bunching" refers to a redistribution or focusing of atoms in an optical potential. For 
the RIR and CARL, these terms are synonymous. If recoil effects are neglected, that is, if the center-of-mass motion is 
treated classically from the outset, the total atomic density is conserved for each velocity subclass of atoms (neglecting 
collisions). A homogeneous atomic density remains homogeneous to all orders in the atom-field interaction. Recoil 



13 



effects allow for a modification of the total atomic density. Whether one calls this "atomic bunching" or "matter 
grating" production is a matter of personal preference f| . The key point is that the modification of the total atomic 
density results entirely from effects related to recoil on the absorption, emission, or scattering of radiation. To 
lowest order in the atom-field interaction, the matter grating or atomic bunching consists of a spatial modulation 
of the atomic density having period 27r/|ki — k2 1 . With increasing field strength, higher order spatial harmonics 
are produced, corresponding to "higher order matter gratings" or "focusing" or "atomic bunching." Of course it is 
possible to derive an effective potential of the form 



(i 



Xi I 



|X2(*)|- 



h 
A 



XiX2(*) e 



i(q-R+<5t) 



XiX2(*)< 



-i(q-R+St) 



ID (ll 



(72) 



without quantization of the center-of-mass motion and to consider classical motion in this potential. However, since 
the effective potential is proportional to h, any changes in the atomic density vanish in the classical limit. This is in 
contrast to bunching in the free electron laser where the effective potential does not vanish in the classical limit. 

Are there situations where matter gratings are produced by fields without any contribution from recoil? The answer 
to this question is "yes," provided one considers the matter gratings associated with individual internal atomic states 
rather than the total atomic density. For example, for an ensemble of stationary, closed, two-level atoms there is a 
spatial modulation in the population difference between the excited and ground states produced by the pump and 
probe fields. Scattering of the pump field from this spatially modulated population difference leads to a dispersion- 
shaped probe absorption profile centered at a probe-pump detuning S = in the limit of large atom-field detuning 
A. The amplitude of the dispersion profile varies as the square of the pump field intensity and its width is equal to 
the excited state decay rate [[L5[ . It is also possible to have gain profiles with a width corresponding to some effective 
ground state decay time if one considers open systems in which spontaneous emission also plays a role. For example, 
if the ground state consists of two hyperfine states and the fields drive transitions between only one of these ground 
states and an excited state, then, as a result of spontaneous emission, both ground state sublevel populations can be 
spatially modulated (although the total atomic density - the sum of all ground and excited state populations - remains 
constant in the absence of recoil). It is possible to monitor the atomic gratings in specific ground state hyperfine 
levels by using radiation that couples only the targeted ground state sublevel to an excited state. 

Finally, I would like to discuss some experimental implications of the RIR and CARL. To observe the spectral 
features of the RIR and CARL, it is best to use cold atoms in a collisionless environment. Moreover, to isolate the 
effects of interest, one must use experiments which involve closed, two-level transitions or signals that depend only 
on total atomic state density. The RIR have already been observed in several experiments involving laser-cooled 
atoms |^-^| . The experiments in which evidence for CARL was claimed [^l],^2| were all carried out under conditions 
(collisions, radiation trapping, large Doppler widths) which do not favor observation of CARL. The results of these 
experiments can be explained by conventional theories in which recoil is neglected 23 . As was mentioned previously, 
the CARL limit is essentially a subrecoil limit. Although the CARL limit has been stated as r = ^ < 1, when Q > 1, 

a better estimate for the CARL regime is ^Q 1 ^ 3 < 1 |lSj ]. For densities of order 10 18 atoms/m 3 , it is possible to 

achieve values of Q 1 ^ 3 of order 100-1000. Thus, it may be possible to observe CARL for atoms cooled to or somewhat 
above the sub-Doppler limit of laser cooling. The linear gain coefficient in this case is of order 10 s , which must 
exceed any cavity loss. 

Is it possible to observe the RIR in a thermal vapor using pump-probe spectroscopy? The linear absorption rate is 
" 22 22. which implies that the ratio of the recoil-induced gain g [Eq. (|47])] to the linear absorption rate is of 



.9i = 
order 



he 
K\ X l\ 



A 2 



where E„ 



. ., — mu /2 is the kinetic energy of the atoms. For sub-Doppler cooled atoms, this ratio can 

be greater than unity, but it is small at room temperature. Still it might be possible to use modulation techniques to 
isolate the RIR contribution to the probe absorption. For the RIR signal to be larger than the background, dispersion- 
like contribution to the probe absorption that varies as |xi I 72/ (A 5 qw), one requires that the ratio -A^ j^-^ 
be greater than unity. This can be achieved at room temperature for sufficiently large A/xi- 

Perhaps the best way to observe CARL would be to use an atomic beam, transversely cooled below the recoil limit. 
The beam can be passed through a cavity with some transit time r . It is not difficult to extend the theory to allow 
for a transit time effects through an effective decay rate T for ground state atoms. For subrecoil cooled atoms, Eq. 
(pq) is replaced by 



(Bit)) 



2Xie 



*„—i5t 



A 



dt' X2 {t') exp(- r +^)(*-*') sin [w,(t - if)] , 



(73) 



Eq. 



by 
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d 2 (B(t)) /dt 2 = -2Td (B(t)) /dt - (oj 2 q + T 2 ) (B(t)) - XlX2( ^ ge , (74) 

and the indicial equation ( pl| ) by 

s 3 + (2r + iS)s 2 + (lu 2 + T 2 + 2iTS)s - i^Q - uj 2 S - T 2 6) = 0. (75) 

With the inclusion of decay, the MR limit, in which (B(t)) depends locally on X2(t), is T ^ Q 1 /^LO q . Thus the CARL 
limit occurs when T < Q x ^ui q . Decay tends to diminish the gain parameter when Q > 1, but can actually reduce 
the gain threshold when Q < 1. For example, if S = 0, the threshold condition is reduced from Q > to Q > if 

rvo. 
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